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Abstract: The behavior of the single-species reaction process A + A — > O is examined near an 
impenetrable boundary, representing the flask containing the reactants. Two types of dynamics 
are considered for the reactants: diffusive and ballistic propagation. It is shown that the effect of 
the boundary is quite different in both cases: diffusion-reaction leads to a density excess, whereas 
ballistic annihilation exhibits a density deficit, and in both cases the effect is not localized at the 
boundary but penetrates into the system. The field-theoretic renormalization group is used to 
obtain the universal properties of the density excess in two dimensions and below for the reaction- 
diffusion system. In one dimension the excess decays with the same exponent as the bulk and is 
found by an exact solution. In two dimensions the excess is marginally less relevant than the bulk 
decay and the density profile is again found exactly for late times from the RG-improved field 
theory. The results obtained for the diffusive case are relevant for Mg 2+ or Cd 2+ doping in the 
TMMC crystal's exciton coalescence process and also imply a surprising result for the dynamic 
magnetization in the critical one-dimensional Ising model with a fixed spin. For the case of ballistic 
reactants, a model is introduced and solved exactly in one dimension. The density-deficit profile 
is obtained, as is the density of left and right moving reactants near the impenetrable boundary. 
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Introduction 



Systems of reacting and coalescing particles, whether chemical or more exotic, are excellent exam- 
ples of statistical systems far from equilibrium. Such processes are widespread throughout nature 
and cover a broad range of physical phenomena. Though the microscopic details of specific reac- 
tion systems may vary enormously, it is often the case that the global behavior of such processes is 
determined by only a few of its properties: the number of particles that interact locally in a single 
reaction and the dynamics of a single reactants' motion. Many of the simple reaction schemes have 
been extensively studied, for example the first-order reaction processes A + A^OoxA + A— > A. 
Also, much is known of the more complicated two-species process A + B — > O where the formation 
of reaction fronts and spontaneous segregation of species Q can occur. These reaction pro- 
cesses have been analyzed for the case of reactants that diffuse and also that have ballistic motion. 
The former is the appropriate choice for systems where the mean-free path is much less than the 
inter-react ant distance, and the latter is appropriate for the contrary case, for example gas-phase 
reactions. Of course, given that the inter-reactant distance is generally decreasing as the reaction 
process continues it can occur that in some physical systems a cross-over will exist from ballistic 
to diffusive behavior |J. 

The majority of studies have considered reaction processes that occur in a translationally 
invariant system, for example in systems with periodic boundary conditions or of infinite extent. 
In this paper we analyze the A + A — > O reaction near a fixed, impenetrable boundary. Such a 
boundary represents the vessel containing the reactants, though as will be discussed, it can also 
represent other phenomena in esoteric reaction systems. 

In section 1 we describe the case of diffusive reactants near a boundary and show that a 
density excess is formed which extends into the system. A concise account of these results has 
been presented elsewhere Q, but without a detailed description of the method. The late-time forms 
of this density excess are found exactly in one and two dimensions and we discuss a surprising 
implication for the dynamic magnetization of the critical one-dimensional Ising model. As a 
counter-point to the diffusive case, in section 2 we consider the case of ballistic annihilation near 
an impenetrable boundary and show that it has the opposite effect: a density deficit is formed at 
the boundary. A model for the ballistic process is introduced and solved exactly with the density 
profiles of left and right moving reactants obtained. Finally, we close with a discussion of the 
results obtained and open questions. 

1 Reaction-diffusion near a boundary 

In this section the effect of an impenetrable boundary on the dynamics of a single-species reaction 
process with diffusive reactants is studied. This universality class comprises the annihilating 
random walk A + A —* O, coalescing random walk A + A —* A and any combination thereof 
||]. Though the results below will be explicitly derived for the annihilating random walk, they 
are identical for every member of the universality class, albeit with a trivially altered pre-factor. 
As well as chemical processes, the dynamics of single-species reactions systems are exhibited in a 
broad range of physical phenomena. For example, the annihilating random walk models the domain 
coarsening of the one-dimensional critical Ising model || . Starting from an initially uncorrelated 
random state at zero temperature, the Ising system evolves by Glauber dynamics. The domains 
increase in size by the process of domain wall annihilation and an isolated domain wall performs a 



random walk (in zero applied field). Such a domain wall can therefore be thought of as a reactant 
in the single-species A + A — > O reaction. (A similar mapping holds between the coarsening of the 
oo-state Potts model and the A + A — > A system). A second example is the exciton coalescence 
reaction A + A — ► A which is seen experimentally in the TMMC crystal |7|, |8[. The non-trivial 
decay exponent predicted theoretically for the laser-induced electronic excitations is seen over 
many orders of magnitude. The boundary conditions to be studied in this section are realized 
in the TMMC crystal with doping by Mg 2+ and Cd 2+ ions that act as perfect reflectors for the 
annihilating excitons |J. 

We first review the relevant results for the reaction process in a translationally invariant un- 
bounded system, which will be referred to as the bulk case. Following that the model will be defined 
and the basic results stated. A detailed description of the method will then be given. The dynamics 
will be written in terms of a master equation and then mapped to a field-theoretic representation. 
The universal properties of the density excess will be obtained and the asymptotically-exact den- 
sity profile found in two dimensions. Motivated by the mapping to the Ising system, a model will 
introduced and solved exactly in one dimension. The section is closed with a summary of results 
obtained. 



Summary of known results for the unbounded system 

A wealth of research exists on the single-species A + A — > O process in unbounded systems 



1 10]— 1 15]. There are many realizations of the dynamics, but the resulting late-time behaviors 
show a strong degree of universality and are independent of the details of particular models. The 
annihilating random walk describes the process whereby diffusing A particles, moving in a d- 
dimensional space, may react pairwise on contact, at a rate A. The principal quantity of interest 
is the decay of the density g of react ants as a function of time t. A starting point in the analysis 
of such systems is the mean-field approach. This corresponds to writing a self-consistent equation 
for the average reactant density ~g, thereby ignoring all spatial correlations 

— = DV 2 p — 2Ao 2 , with the late time result g ~ — — (1) 
at ' ' 2 Xt 

where we have introduced the diffusion constant D. Neglecting the effect of correlations is equiva- 
lent to assuming that the reactants remain well mixed throughout the reaction process. However, 
due to the statistics of random walks in two dimensions and below, simple diffusion of particles 
itself is not sufficiently fast to maintain a well-mixed state. Reactants that are close together react 
fast leaving reactants that are more widely separated, i.e. anti-correlated in space. For this rea- 
son, the mean-field result (jlj) loses its validity below two dimensions. Also, in these lower spatial 
dimensions the fact that random walks are recurrent means that reactants come into contact many 
times, each time providing an opportunity for a reaction to take place. This implies that even a 
small reaction rate A does not limit the global rate of reaction. In two dimensions and below the 
reaction becomes diffusion limited and is independent of the parameter A. The correct forms for 
the density in two dimensions and below are given in Table [l] at the end of this section. 



Reactions near a boundary - new results 

We now define the model that will be analyzed throughout the remainder of this section. The 
dynamics of the A particles are the same as for the unbounded case except that now the geometry 
is restricted by the imposition of a hyperplane boundary, confining all particles to the positive 
half-volume. To be specific, the model is defined on a hypercubic, d-dimensional lattice which for 
convenience has a lattice spacing of unity. This (i-dimensional lattice is infinite in d — 1 transverse 



2 



dimensions and semi- infinite, sites 1, 2, • • • , oo in what will be called the z direction. Initially, the 
lattice is filled with an average density qq of A particles that have two components to their behavior, 
diffusive motion and mutual annihilation. They can diffuse randomly throughout the lattice by 
hopping at a rate D to any neighboring site, as long as the restriction z > 1 is maintained. The 
diffusion is independent for each particle and hence multiple occupancy of a lattice site is allowed, 
leading to bosonic statistics. However, if there are n > 2 particles on a site, a reaction can occur 
there with a rate Xn(n — 1) reducing n by 2. 

In the same way as for the bulk case, the dynamics of the model can be approximated by a 
mean-field description of the same form as equation (|l]) but with the added restriction of zero 
current at the boundary. This amounts to enforcing a vanishing density-gradient at z = which is 
a restriction compatible with the unbounded solution. Hence, the mean-field equation asserts that 
the boundary has no effect. Nevertheless, as was shown above mean-field cannot be relied upon 
in two dimensions and below because of the importance of correlations. In fact, even a simple 
argument shows that a density excess develops near the boundary in low dimensions. Consider 
the dynamics of the model in one dimension, up to a time t and far from the influence of the 
wall. Since random walks are recurrent in these dimensions, in the bulk most particles within a 
diffusion length Ij, ~ \/2Dt will have interacted and annihilated. This leads to a density in the 
bulk of the system of Qb — 1^ = c\,j\ft. Close to the wall, within a distance of the order of y/Dt, 
the diffusion length is smaller. The density near the wall is then g w ~ Z" 1 = c w /\/t, with c w >Cb- 
As 

Cw — e w {z / Df) the argument implies that there is a density excess near the boundary which 
propagates into the system diffusively. 

This is a qualitative argument, however we will show analytically that in two dimensions and 
below this is indeed the case: the reactant density has the form of a constant background density 
qb given by the well-known unbounded results in Table [I], and a fluctuation-induced density excess 
QE 

g(z,t) = Qs(t) + qe(z, t). (2) 
As will be demonstrated, the density excess has the universal dimension-dependent form 

qeM = ^d^ fd {-im)- (3) 

In the following two sub-sections we use both the field-theoretic renormalization group and an 
exact solution to determine the asymptotic forms of the scaling functions fi and f\ for two and one 
dimensions respectively. A reasonable amount of technical details are given, but all the important 
results derived can be found in the summary at the end of this section. 



1.1 Two dimensions and below 



In the following we will examine the late-time scaling behavior of the lattice model described above. 
First, the dynamics will be written as a bosonic master equation and mapped to a continuum field 
theory via the coherent-state formalism. The theory is then regularized and the renormalization 
group applied to obtain the universal properties and non-perturbative results valid for late time in 
dimensions two and below. The method is standard Jl5|, [l6|, 17], but is technically more compli- 
cated than the unbounded case 



14 1 due to the lack of translational invariance that necessitates 



calculations in real space. 



The field theory 
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The first step in the field-theoretic calculation is to write the dynamics of the model in terms of a 
master equation. This describes the flow of probability between microstates of the system defined 
by the set of occupation variables {rii} where is the number of particles on site i. Due to the 
bosonic statistics of the A particles it is convenient to write the master equation in a second- 
quantized form. Each configuration, {n{\ is assigned a vector in a bosonic Fock space | with 
|0) denoting the no-particle state. Using and a k as the bosonic operators for site k the state 
ket of the system \P(t)) is written 

\P(t))=J2 P ({n j }:t)Y[(alr\0). 

M * 

Here P({rij} : t) is the probability of being in a particular microstate {fij} at a time t, the sum is 
over all microstates, and the product is over all lattice sites. The master equation for the model 
can now be written in the form dt\P>= —TC\P> where the evolution operator Ti is given by 



« = £ 



D^4( ai - % )-A(l-(a;) 2 )a 2 



The sum i is over all lattice sites, and the sum j is over all of site i's neighbors, with the condition 
that both sums are restricted to the half-space. We now introduce two objects, the projection 
state (| = (\a* and the initial state \qq) with average particle density go, that allow observables 
such as the density Qk(t) to be written 

(| = (0|e^ and \g ) = e* ^^" 1 ^}, 
6k(t) = (\a k e~ nt \g ), 

where the sums are again over all lattice sites j with z > 0. The field-theoretic representation 
is obtained using the coherent-state formalism [jll| . The continuum of microstates of the system 
are now given by configurations of the field variables {4>q,4> } which are analogous to {a, a) — 1}. 
Observables are written as integrations over the appropriate weighting function <So(<fo> ^o) which 
we will call the action. For example, the density at position r and time t is written 

g(r,t) = J Vct> V^ Mr,t) e- So ^\ (4) 
5 o (0o,0o) = S D + 2\fdt( al d r^4>l + X f alt f al d rj\<& - g f d d rfo(0). (5) 

JO Jz>0 JO Jz>0 Jz>0 

The diffusion component of the action So defines the propagator 



g(r f ,n:t) = [Gizf-Zi-.q + Gizf+Zi-.tyHGiyf-Viit), (6) 

M 

where we have used the notation r = {z,yx,y2 ■ ■ •}, with the set {y} representing the dimensions 
parallel to the boundary and the function G(x : t) is the Gaussian 

1 ( x 2 \ 

G(x : t) = —-pr exp 



(A^Dt) 1 ' 2 "\ Wt t 
Regularization of the theory 
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(a) 



2gD[f 



2(Z g -l)gD\l £ 



(z g -i)gD\i 



(b) 




y 

+ — = finite 



Figure 1: The diagrammatic representations of (a) the interactions and counter terms (with blob) 
and (b) the one-loop regularization. The propagator is denoted by a straight line with an arrow 
signaling the time direction. In the diagramatic equation (b) Pi and P c j are the one-loop and Zg 1 ^ 
counter-term contributions respectively. 

From power counting, it is seen that the upper-critical dimension d c = 2 of the theory is unchanged 
from the bulk case. The theory contains divergences, by virtue of the continuum limit and to 
write the density in a calculationally useful form we need to replace all the bare quantities with 
renormalized ones and render the theory finite. This is done using dimensional regularization in 
d = 2 — e. In fact the renormalization of the theory is straightforward as there is no propagator, 
initial-condition or field renormalization: only the reaction rate is renormalized. Because the 
propagator is not dressed by the interactions, Eq. @ is the full result showing that the boundary 
remains effectively reflecting on all scales. In the language of surface critical phenomena, this 
corresponds to the special transition Jl9| persisting at all orders and is different from the behavior 



frequently seen in equilibrium surface critical phenomena [20, 21] and in related non-equilibrium 
systems f|^, 23]. 

The bare action So = S + S c t is now rewritten in terms of the renormalized action S and the 
counter-term action S c t- We also introduce the dimensionless interaction parameter g defined by 
A = Z g gDfi e with [i arbitrary but with dimensions of inverse length. Thus 



S D + gD^i 
(Z g - l)gDfi 



dt 



d d rdxj> 2 + 



z>0 



dt 



z>0 



d d r$$ 



00 



2>0 



d d r(f> 



(7) 



2 / dt 

o 



d d r66 2 + 



2>0 



dt 



z>0 



d d r4?(i) 2 



The diagrammatic representation of the two reaction terms and the counter terms are given in 
Fig. 0(a). 

can now be calculated up to the one-loop 



The renormalization constant Z n 



1 + zP + 



order. The minimal-subtraction regularization is done perturbatively (in the absence of the initial 
conditions) by demanding that calculated quantities at a given order of perturbation theory are 
finite. As an example, consider the correlation function 



P = {<f>(c, t)<f>(a, 0)0(6, 0)) = [Pi + P a ] + 



finite 
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where (...) denotes an average of the fields over the action (^) without the initial-condition term. 
We have used the notation P\ to denote the contribution to this correlation function at the one-loop 
level and P ct is the appropriate counter-term, see Fig. 0(b): 



P l = 8gD 2 /j 2e f dh f d d n f dt 2 f d d r 2 

JO Jz>0 JO Jz>0 

G(c-ri,t- t 1 )g 2 {r 1 - r 2 , h - t 2 )Q(r 2 - a, t 2 )Q{r 2 - b, t 2 ), 

Pa = -Z^HgDpf f dh f d d na(c-ri,t-ti)g(ri-a,ti)g(ri-6,ti). 

Jo Jz>0 

Using the above equations the scheme requires that [Z^ — 2Dg/j, e K] is finite, where K is defined 
through 

h * f jd„ 0(r 2 -a,t 2 )G(r 2 -b,t 2 ) 2 



Jo Jz>o G[ri - a,*i)^(ri - b,tx) 
The singular part of K can be extracted so that 

2Dg^K = Z^ = JL + 0(e°). 

Hence using the renormalization constant Z g = 1 + gj^Lixe) the one-loop beta function is 

0(9) = |-( 5 -2vre) ind=2-e. (8) 



The fixed point structure is unchanged from the bulk result found in [14]. This is understandable 
as physically the renormalization of A is connected to the fact that random walks are recurrent in 
two dimensions and below: a feature unaffected by the presence of a boundary. With the theory 
regularized, it is now possible to proceed with a perturbative calculation of the density that is 
valid for early times 

The perturbative density 

To calculate the density g at early times, three objects are needed: the tree-level contribution g(°\ 
the one-loop contribution and the counter-term 

g = e (°)+ e « + e g>. (9) 

All these are calculated using the renormalized action with the initial conditions included 

g(r,t) = ( V<pV~4> <j)(r,t) e~ s ~ Sct . 



where S and S c t are defined in Eq. (R). However, as will be justified later, an expansion is made 
in the parameter g$ 1 with only the leading term retained. This is equivalent to the limit of high 
initial density. 

(i) The tree-level contribution g^ 

The tree-level density is obtained by the sum of all diagrams given in Fig. ^(a) with the sum 
rewritten to yield an integral equation. The result is found to be equivalent to the bulk case, as is 
expected from the mean-field analysis earlier. 



Q 



(0) 



Q0 1 



1 + 2g g^Dt 2gfi 6 Dt 



+ O(g 1 ). (10) 
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(0 

p - - 

1 1-loop 

Figure 2: The diagrammatic representations of (a) the tree-level density calculation, (b) the clas- 
sical propagator and (c) the perturbative density to the level of one loop, which reading from left 
to right are g(°\ g^ and g^) 




Before proceeding to the one-loop calculation, the dressed, tree- level propagator must be calculated. 
Following [|l4| we call this the classical propagator Qc and it is defined through 



Qc = (<Hr-i,ii)<Hr 2 ,t 2 )> 



c 



with only the tree diagrams included in the average, see Fig. |2|(b). This propagator is evaluated 
in real space to give 



(ii) The one-loop contribution g^> 

Using the classical propagator the one- loop contribution takes a simple form, see Fig. §(c), 

gW = 4DVV / dfinPdhf d d r 2 gc(r-r 1 ,t-t 1 )gUri-r 2 ,t 1 -t 2 )(g^(t 2 )) 2 . 

JO Jz>0 JO Jz>0 K ' 

To simplify the notation we introduce Gb{t,£), the bulk propagator, defined through 

GB(r f ,n:t) = Gizf-Zi-.QUGiyf-Vi-.t). (12) 

M 

In the large go limit the one-loop contribution then takes the form 
gW = _L J* dtl J* 1 dt2 (^j (J B + Ie ) } W ith 
/b = / d d ri J d d r 2 GB{r-r 1 ,t-t 1 )g 2 B {r l - r 2 ,h - t 2 ) 
Ie = J dfinj d d r 2 g B {r-r 1 ,t-t 1 )g B {ri-r 2 ,t 1 -t2)9B{ri-r' 2 ,tx-t2)- 
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Note that the half-volume integrations have been exchanged by integrations over the volume. 
The integral has been separated into two parts Ib and Ie corresponding to the bulk and bound- 
ary/excess contributions. The notation r' 2 denotes a set of coordinates in which the z component 
is minus that of r 2 , but all other coordinates are unchanged. We consider these two quantities Ib 
and Ie separately. 

Firstly, the bulk integral Ib is performed and gives the following one-loop contribution to the 
perturbative density 



Ib 



1 



2(27Tl?(ti - t 2 ))V5 



= (^(ik + ii( 2l0g ^- 5 )) +0( '»- 



This result is the same as that found in |14] and provides a uniform background density. 

Secondly the excess contribution Ie is considered. The integrations over the d — 1 dimensions 
parallel to the surface yield the result as the integral Ib (with the power d replaced by d — 1). 
However, the integration over the direction perpendicular to the wall yields new a z dependent 
result. Combining these results one obtains for Ie 

1 1 f 7 2 \ 

J E = 77 r^,,/o T, , Ul-^HTC, \ TTTTt ex P 



{8irD) d / 2 (t 1 -t 2 )( d - 1 )/ 2 {2t-t 1 -t 2 ) 1 / 2 \ 2D{2t-t 1 -t 2 ) 
Using this we obtain the contribution to g^\ 



ri) 1 ( 1 \ . / z 



2 



°* = <j5WA^) h Ul* +0 (e)) where 



z 2 



/2 '^J = Jo ds L dq {s) ( s - q y/ 2 (2- s - q y/ 2 - (13) 



We were unable to evaluate this double integral, except at z 2 /Dt = where the result /2(C)) = 
3/2 + 7r — 37r 2 /8 was found. However, it has the asymptotic form f 2 {^ 2 ) ~ exp(— £ 2 /2)/£ 3 which 
is valid for £ 3> 1. 

(iii) T/ie counter term 

As before in terms of the classical propagator this contribution takes the simple form, see Fig. 
1(c). 

rt 

J, I Kir-: , , , I ..MM, f | 



JO Jz>0 v 

In the large go limit this is evaluated in a manner similar to g^ to give g^ = —Zg 1 ^ j2Dtg\f . 



All the required perturbative results have now been obtained up to the order of one loop. 
These tree-level and one-loop results are combined and rewritten as in equation @ as the sum of 
a spatially independent, background density gs and a spatially dependent, excess density qe 



QB{t) 



QE(t) 




1 



4areDtfJL ( 



+ G(g) 



87T(Dt) d / 2 



(14) 



(15) 
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where the results here are correct at the level of one-loop. In the next section this perturbative 
density will be used to obtain the late-time results by the use of the renormalization group. 



The flow equations 

A Callan-Symanzik equation for the density can be obtained by noting that any observable must 
be independent from the arbitrary /i. Using this fact and dimensional analysis allows one to write 
two partial differential equations which can be combined to give 



d d d . d 

o log z a log Dt o log go og 



g(z,Dt,go,g,fi) = 



(16) 



d 



+ d 



0.. 



d log s 

where we have rewritten the sum of partial derivatives as a total derivative with respect to a 



scaling variable s. Following the notation in [14] we write X = X(s) for early-time quantities and 
X = X(l) for late-time quantities, 



s g 



z = sz, t = s t, g = s g , g 



(3(g)dlog(s), 



(17) 



and obtain the exact relation between a density with arguments {z,Dt, ^cbflSA*} to a density with 

— \ g(z,Dt,g ,g,(j,). 



g(z,Dt,g ,g,fj) 



(18) 



Next we will take the limit of large t/t (s — ► 0). In this case the equation relates a system at late 
times t to the density of a system with a large initial density go measured at early times t. 



The late-time density excess below two dimensions 

The scaling equation derived above is now used to obtain the late-time density. Using the notation 
for early-time quantities defined above, the perturbative excess, Eq. (15) is rewritten with t — * t 
etc and inserted into the scaling equation (jl^)- In dimensions d = 2 — e the quantity g tends to 
its fixed point value g = 2ire as the ratio t/t increases. The spatially independent and dependent 
density components become in the late-time limit 



QB(t) 
QE(z,t) 



1 



(Dt) d / 2 
1 

8TT(Dt) d / 2 



1 | 21og(8^) 
Aire 167r 



h 



z 



2Dt 



+ 0(e). 



(19) 
(20) 



As expected, the bulk value is unchanged from the well known case derived in [14] and given in 
Table |l[ However, also present is an excess gE that decays at the same rate as the bulk case, i.e. 
it has a fixed amplitude with respect to gs- It shares the same universality as the bulk case in 
as much as it is independent of £0 and the reaction rate A. It should be noted that the excess is 
not localized at the boundary but extends diffusively into the system by virtue of the functional 
dependence of on the ratio z 2 /Dt. 



The late-time density excess in two dimensions 

In two dimensions the beta function takes a different form (3 g cx g 2 . In this case we have the scaling 
of the reaction parameter in the perturbative density g ~ [4-/r log(t/t)]" 1 . The leading order terms 
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in time are now taken, which involves the RG-improved tree-level calculation for the bulk density 
and the one-loop calculation for the excess density 

latfDt/Efy (1\ log(i) (I 
eB(t) = 8,Dt + °{- t ) = 8^ + °U 

= ^Dt h (^) + °{TTk))- (21) 

Again, the bulk density is the known result given in Table |l|. Now the leading-order term of the 
excess is the fully universal, asymptotically exact result for two dimensions. We will discuss this 
further in the summary. 

1.2 Exact solution in one-dimension 

The renormalization group has provided the universal properties of the density excess as a function 
of dimension and the late-time form was found exactly for the case of two dimensions. However, 
the epsilon expansion is not expected to give accurate results in one dimension. In this section we 
solve the one-dimensional case with an infinite reaction rate. However, by virtue of the universality 
demonstrated above, the results to be described are also valid for any finite reaction rate. 

The model is now defined in dimension d = 1 with an infinite on-site reaction rate. Since the 
reaction rate is infinite, each lattice site can be occupied by at most one particle. Denoting a 
particle on site k by A\. and an empty site by Ok the dynamics of the model are then 

OkAk+i <-> AkOk+i with a rate D 

A k A k+ i -> OfcOfc+i with a rate 2D, (22) 

where k is restricted to the positive integers. At t = the state of the system is chosen to be an 
uncorrelated random initial state with density 1/2. 

Homogeneous reaction systems in one dimension have been solved exactly by various methods 



10]-[O]. For inhomogeneous systems, a general method was recently proposed in which the model 



is solved by mapping it onto a dual system ]2j, |25|. Using this method, the master equation is 



first written in a quantum spin chain representation. Each configuration is represented by a vector 
|si, S2, . . .) = where Sj takes the value 1/2 if site i is empty and —1/2 is site i is occupied. 

The state ket of the system is defined similarly to the bosonic case to be 

\P(t)} = ]TP({s;} :t)\{ 8i }} d t \P> = -H\P> 
W 

oo 

n = ~ D E (4%fi + *Ui*k + 2 44+i - a - **)/2) • (23) 

fc=l 

Here = (of ± iot)/2, where the cr's are Pauli matrices. Averages are calculated in the same 
way as the bosonic case, by using the projection state (| = Z){s}({ s i}l- It was shown that using 
a similarity transformation the system can be mapped onto a dual process with the following 
dynamics 

AkAk+i — ► OfcOfc+i rate 2D for sites k = 1 • • • oo 

OkAk+i — ► AkOk+i rate D for sites k = • • • oo (24) 

AkOk+i — ► OkAk+i rate D for sites k = 1 ■ ■ ■ oo 

^1^2 -> 0i0 2 rate D. 
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Denoting the dual Hamiltonian as TC, the required density takes the form 



g k (t) = -{2\e- nt \k-l,k). 

Here (2| is the sum of all two-particles states and \k— 1, k) is the initial state with only two particles 
located at sites k and k — 1. However, since the reaction rate is infinite the two-particle transition 
probability can be written in terms of single-particle transition probabilities, giving the density as 



oo oo 

(*) = o T, E [T m ,k-l(t)T n>k (t) - T n ^ x (t)T m>k (t)] , 



(25) 



m=0 n=m+l 



where T n ^(t) is the probability of a particle starting at site n to be at site k at time t in the dual 
system with the forms 

T z , a {t) = B{z-a:t) + B{z + a:t) 

oo 

m=l 

Here B(m : t) = e~ 2Dt I m (2Dt) is the single-particle propagator with I m the order-m modified 
Bessel function. Using these, results the double sum (|25|) is reduced to 



Qk(t) 



1 

2 L 



B(l -k:t)J2 B ( m ~ k:t )+J2 B ( m ~ k:t ) i B ( m ~k:t) + B(m + 1 - k : t)] 

m=3 m=2 

oo oo 

B{k : t) B (™ + fc : t) - £ B(m + k : t) [B(m + k : t) + B(m + l + k:t)] 

m=2 m=l 

oo oo 

B(k : t) B ( m -k:t)+B(l-k:t)J2 B ( m + k : t) + B(k : t) + B(l — k : t) . 



m=2 



m=l 



Taking the late-time limit, which is valid when V Dt is much greater than the lattice spacing and 
writing k = z, we obtain 



g(z,t) 
g(z,t) 



1 



V87rDt 
1 



erf {^ J^ D ^ j + v^erfc ^ 
+ QE{z,t), 



k 



exp 



ADt 



V8irDt 

where again the density is split into bulk and excess qe components 



QE(z,t) 



1 



/i(£ 2 ) = >/2erfc(-i= 



y/torDi \2Dt) 
exp 



(26) 



(27) 



erfc (£ ) . 



As expected the bulk component of the density qb is found to be identical to the infinite, unbounded 
case [12, [DJ. The asymptotic form of the function is /i(£ 2 ) ~ ex P( — £ 2 /2)/£. 
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d 


late-time bulk density qb 


late-time excess density qe 


1 


i 

(87r£»t) 1 /2 


(87rDt)V2 ( 




) 


2 - e 


1 


[1 + | (21og(87r)-5); 


+ 0(e) 






)+0(e) 


47re(Dt) d / 2 


(87rDt) d / 2 ^ 2 


2 


log(t) 


(87rDt)^ 2 






> 2 


1 

2At 


sub-leading 



Table 1: The known results for the late-time bulk density and the new results for the bounded 
system. 



1.3 Summary of the diffusive case 

Results have been presented from an analysis of a reaction-diffusion process near an impenetrable 
boundary. It was found that in two dimensions and below a density excess forms near the boundary 
and extends into the bulk diffusively. The density can be written in the following form 

g(z,t) = g B (t) + QE(z,t), 

where QB{t) is the bulk result seen in a translationally invariant system and qe is the boundary- 
induced excess density - see Table |l[ Using the field-theoretic renormalization group the universal 
properties of the density excess were found. It was shown that the excess is independent of both 
the initial density and the reaction rate A. In two dimensions it was found that the density excess 
is marginally less dominant than the bulk, with the following asymptotic (large z 2 /2Dt) limit 

1 /2Dt\ 3 / 2 ( z 2 \ 
QE = ^ m {—) exp [- m m two dimensions. (28) 



The full asymptotically-exact late-time result is given by ( |2T| ) with the double integral ( |13D and is 
plotted in Fig. 0. 



In one dimensions the excess density was obtained exactly (27) by the solution of a corre- 
sponding quantum spin system and is also plotted in Fig. [||, with the asymptotic (large z 2 /2Dt) 
limit 

1 /2Dt\ 1/2 ( z 2 \ 
QE = ^[-im) -one dimension. (29) 

Here we see that the excess and bulk density share the same leading order decay of g ~ t~ l l 2 . The 
ratio between the density at the boundary and the density far in to the bulk below two dimensions 
is thus independent of all system parameters except the dimension of space. In one dimension it 
is numerically equal to v2 and it is interesting to note that the renormalization group gives quite 
a fair indication of this ratio 

which is ~ 1.47 for one dimension (e = 1). It should also be noted that above the critical dimension 
there is a sub-dominant density excess. However, these effects are transient and quickly decay to 
yield the predicted mean-field result. 
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Figure 3: The density excess in one and two dimensions, f\ and fi- The renormalization group 
results are compared with simulations at time t = 10 (o) and t = 80 (*). The exact results /i 
calculated for infinite reaction rate is compared with simulations for finite reaction rate A = 1/2 
(+) demonstrating the universality. 

We end the discussion of the diffusive case by considering the implications of the solution in 
d = 1 to the coarsening dynamics of the critical Ising model. As explained earlier, the domain wall 
dynamics in this Ising system can be mapped onto the mutually annihilating random walk ||. The 
boundary in the reaction system corresponds to a fixed spin in the magnetic system. Consider the 
magnitude of the coarse-grained magnetization near this fixed spin. This quantity is a function of 
the local density of domain walls - the fewer the domain walls the higher the magnetization. The 
density excess (|27]) implies a higher density of domain walls near the fixed spin. This yields the 
counter-intuitive result that the absolute value of the coarse-grained magnetization is lower near 
fixed spin. 
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2 Ballistic annihilation near a boundary in one dimension 



In the previous section the effect of an impenetrable boundary on the single-particle reaction- 
diffusion process A + A —* was examined. It was found that for diffusing reactants a density 
excess was produced at the boundary In this section we will demonstrate that this effect is not 
independent of the reactants' dynamics: for the case of ballistic annihilation a density deficit is 
formed. 

By ballistic motion, it is meant that reactants travel in straight lines until they collide and a 
reaction occurs. These dynamics are relevant for gas-phase and some excitonic reactions j| and 
a good review can be found in ||. Because of the deterministic dynamics of individual reactants, 
the only sources of noise are in the initial distribution of particles or if the reaction rate is finite. 
Though a great deal of analytic work exists on these systems, they have largely examined two kinds 
of initial conditions in one dimension. One form of initial condition was introduced [26] to examine 
the dynamics of reaction fronts and has been demonstrated to be universal with respect to a finite 
reaction rate pTj] . However, the case that will interest us in this section is that of random initial 
conditions whereby particles are distributed upon a line with each particles' velocity drawn from 
the same distribution. A broad variety of initial- velocity distributions have been studied (8|, 28, 29] 
and universal effects analyzed for the continuous case However, in these cases the systems 

that have been examined have been infinite and translationally invariant. 

In this section we will generalize the binary-velocity model introduced by Elskens and Frisch 
[31]. The behavior of this translationally invariant system will be recalled below and the general- 
ization to the case of a semi-infinite system bounded by a reflecting wall will be introduced. The 
basic results will be stated and then details of the method given. Finally, the section is closed with 
a summary of the behavior in the late-time limit. 



Summary of known results for the unbounded system 

The model introduced in [31] comprises an infinite line upon which particles move with fixed 
velocity +c or — c. When two particles meet a mutual annihilation always occurs, and both 
particles are removed from the system. The initial conditions provide the only source of noise and 
are defined such that at time t = particles are equally spaced along the line and given either of 
the two fixed velocities with equal probability. A mean-field approach predicts that the density 
decay should vary as g ~ t . Nevertheless, the density decay derived from a full solution of the 
model has the form g ~ t -1 / 2 . This is the same exponent as for the diffusive case, but it should 
be stressed that whereas in the diffusive case the changed exponent comes from a dynamic effect, 
in ballistic annihilation it is a property only of the initial conditions. This can be seen from the 
following qualitative argument. Consider the density fluctuations in a domain of length t. The 
typical fluctuation in the difference of particle numbers with v = +c and v = —c inside this domain 
is oc i 1 / 2 . After a time oc £ only the residual fluctuation will remain, so that the density g will 
be oc I- 1 / 2 . Re -expressing £ as a function of t one obtains the correct density-decay exponent 
gcxt-V 2 . 



Annihilations near a boundary - new results 

A reflecting boundary is now introduced into the model defined above. The boundary is placed at 
position z = 1/2 on the one-dimensional line and particles are placed, equally spaced at positions 
z = 1,2, ■ • -. When referring to particle k it is meant that this is the particle that was at position 
z = k at time t = 0. Each particle is given a velocity of either +c or — c with equal probability. 
The velocity of each particle is fixed unless it reaches the impenetrable, reflecting boundary. In 
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this case a particle with velocity — c rebounds with a fixed velocity +c. Reactions occur with 
a probability of unity on contact, so that once the trajectories of two particles meet a mutual 
annihilation always occurs, see Fig. [5] for an illustration of the dynamics. 

The quantity of interest is the density g averaged over all possible initial velocity distributions. 
It will be demonstrated below that this density can be written as a sum of a bulk, background 
density qb equivalent to that calculated in [31] and a boundary-induced density deficit qd 



g(z,t) = g B (t) + QD(z,t). 
The scaling form of the density deficit, given exactly in the following sub-sections, is 

g D (z,t) = — —r-, g (4) » ( 31 ) 



(ct) 1 / 2 * \ct 

and should be compared to the diffusive form which decays with the same exponent ~ t~~ 1 / 2 but is 
a function of the dimensionless quantity z 2 /Dt. In the ballistic case it is also seen that the effect 
penetrates from the boundary into the bulk, though here it is ballistic by virtue of the dependence 
of g on the dimensionless form z/ct. 

A detailed description of the exact solution now follows. First, the probability of any specific 
pair of particles to react is calculated. The method of solution is combinatoric and related to 
the well-known ballot problem |3^] of probability theory. From this fundamental pair-reaction 
probability all other quantities are obtained, including the density profile and the profiles for left 
and right moving particles. 

2.1 Exact solution in one-dimension 
The pair annihilation probability 

Because the model defined above has a reaction probability of one, only mutual annihilations 
between even and odd particles are possible. There are two cases that need to be considered: 
either the even or the odd particle can be nearest the boundary. We obtain all the possible pair- 
annihilation probabilities by considering an odd numbered particle k annihilating either with even 
particle k + (2n + 1) or even particle k — (2n + 1) where n = 0, 1, • • •. 

Consider first the probability of the odd numbered particle k to mutually annihilate with 
particle k + (2n + 1). This probability has two contributions: (i) direct annihilation with particle 
k having initial velocity v = +c and (ii) indirect annihilation with particle k having an initial 
velocity v = —c but rebounding from the wall. For direct annihilation the calculation is identical 
to the homogeneous case [pi]] . A mapping exists to the problem of a random walker returning to 
the origin after 2n + 2 steps for the first time - this is the so-called ballot problem. The probability 
for this to happen, given that the initial velocity of particle k is +c is 



2n 



Next, consider the contribution from the indirect case whereby a particle is first reflected from 
the boundary before annihilating. In this case the initial velocity of particle k is v = —c with 
the total probability given by the product of Eq. (^|) and the probability that particle k reaches 
the boundary. The latter can be found by noting that particles labeled by an even index with an 
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k k+1 k+2 k+3 k+4 k+5 



k+(2n+l) 



©000©© © 

1 1 1 1 1 1 1 

k k+1 k+2 k+3 k+4 k+5 k+(2n+l) 

Figure 4: An illustration of the mapping to the ballot problem. The arrows, representing the 
direction a particle is moving, are considered as moves of a random walker as illustrated above. 



initial velocity v = —c are always annihilated before arriving at the wall. It is therefore equivalent 
to the probability of a random walker not to have returned to the origin after k — 1 steps, 

/l\*-a 5ri x fk-l\ /l\k-2/fc-2\ , , 

where this is conditional on particle k having an initial velocity v = — c. Combining the two results 
of Eqs. ( |32"| ) and (^) gives the net probability that particle k interacts with particle k + (2n + 1), 
Pk,k+(2n+l) regardless of particle ks initial velocity 

1 /l\2n+2/2n + l\ / /l\k-2fk-2\\ , 

W«i = sr+iCs) („ + ij( 1 + y (^)J- (34) 

The final case of the odd numbered particle k to interact with particle k — (2n + 1), Pk,k-(2n+l)i 
can also be derived using the arguments given above 



1 /l\2n+l 2n + 1 . 
V^D-^ls) n + l'- (35) 



The Eqs. (|34j ) and ( |35| ) give all the possible pair annihilation probabilities in the system, where 
it should be remembered that particle k is odd numbered. 



The density calculation 

To ease the notation in this section we set c = 1, corresponding to a trivial redefinition of time. The 
initial conditions of equally spaced particles introduces a discretization of time into the problem. 
We consider the model at integer time steps when particles will be found only at positions z = 
1, 2 • • • which we call lattice sites. In the calculation of the density two cases must be considered: 
(i) z > t when the presence of the boundary does not affect the density profile, (ii) z < t where 
the boundary does affect the profile. 
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Figure 5: An example of the initial conditions studied. The numbers mark lattice sites, and dashed 
lines represent the trajectories of the particles. 



(i) The case z >t 

Here the density is the same as for the infinite-system case |pj| . The contribution to g(z,t) can 
arise from particle z + t and z — t. Using the mapping to the ballot problem we find 




which is the probability of particle z—t not to interact with particles up to and including z+ (2i — 1) 
added to an equal contribution from particle z + 1. Note that the results are valid for t = z since 
the contribution from a particle initially at site 1 with velocity v = — 1 which rebounded off the 
wall is equivalent to the contribution from a "ghost" particle initially at site with velocity +1. 



(ii) The case z > t 

There are two contributions to the density at site z, one from particle z + t with an initial velocity 
v = — 1 and another from particle t — z + 1 with an initial velocity v = — 1 which rebounded from 
the wall at time t — z + 1/2. First, the contribution from particle z + 1 is considered with z odd. 
This particle reaches site z if it does not annihilate with any of the z + t — 1 particles to its left. 
When z + 1 is odd (t is even) we obtain for this probability 



l\*+t-i 
2 



(i) 



z+t-l 

E 

x=2,4,... 



Z + t-l 



z + t-r 

z+t-l+x 
y 2 / 



(5) 



z+t-l 



z + t-T 

z+t-3 



When z + 1 is even (t is odd) the particle never reaches the wall so that the probability of a particle 
initially at site z + t of surviving up to site z is the sum of all the chances of being annihilated 
after site z 
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Figure 6: An example of the exact (a) density profile and (b) densities of left and right moving 
particles at time t = 30. The difference between odd and even sites becomes increasingly less 
pronounced at later times. 



Now, the contribution from particle t — z + 1 is considered also with z odd. This is non-zero only 
if t — z + 1 is odd (t is odd), and is then given by, 



/l\t- z (t- z- 1\ /l\2«-l if, x ( 2: 
UJ 1 4=1=2 )\2) x J^ Yz 



2z+x 



,i\t+z-i(t-z-i\(2z-i\ 

- y )[z-i)- 

Combining the above results from particle z + 1 and particle t — z + 1 we obtain the density at odd 
sites z and at even and odd times. Using similar argumentation allows the results for even sites 
also to be found. The forms for the densities are 



g(z dd, ''even ) = e( todd) = ( o ) 



l\z+t-i{z + t - 2\ 

z+t-3 ' 



/l\z+t-l 

Q(z dd, todd) = Q{ Zeven ; ^etien ) 



z + t-l/z + t-2\ t - z + 1 A - z\ (2z - 2 N 



t-z- 1\ /2af — 1 
^ U-l. 



These results taken together represent the complete density profile at all times. An example is 
shown graphically in Fig. [6|. It is clear from these results that a density deficit is indeed formed 
at the boundary. It should also be commented upon that the even/odd effect, though strong at 
early times, becomes increasingly less significant at later times. We now consider how this density 
is distributed between left and right moving particles. 



The density of left and right moving particles 
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For the case of the unbounded system the density of left and right moving particles is equal by 
symmetry. For this reason it is only necessary to consider the case z < t because the boundary 
has no influence for z > t. Using similar method to the above, the following forms are found 

The density of left-moving particles 



£ , {,Z<>dd; teven j 

£{z dd, todd) 
/l\*+t-l 

V2/ 





(37) 



The density of right-moving particles 

T^\^oddj t e ven) 

7t(z odd , t dd) 



As an example, the left and right densities for the even time of t = 30 are plotted in the second 
graph of Fig. [f| From this graph it can be clearly seen why the deficit occurs. As time progresses 
the remaining particles tend to congregate in groups which move in the same direction. Once such 
a group hits the boundary it mostly annihilates within itself, reducing the density of returning 
(right-moving) particles. 

2.2 Summary of ballistic case 

The effect of an impenetrable, reflecting boundary has been examined on a one-dimensional realiza- 
tion of ballistic annihilation. The system evolves deterministically from an initial random velocity 
distribution and various quantities, averaged over different realizations of this velocity distribution, 
were calculated exactly. In particular the density profile of the reactants as a function of time and 
distance from the boundary was obtained. It was found that, in contrast to the diffusive case 
studied in the previous section, a density deficit gr>(z,t) was formed at the boundary that extends 
into the system ballistically 

g(z,t) = QB(t) + Qd(z, t) 

In the late-time limit, the scaling function g takes the form given in table |2| and the late-time 
form of the full density is shown in Fig. [j]. Also calculated were the density profiles for the 
left and right moving particles, of which the asymptotic functional forms are given in table | 
and are plotted in Fig. ||[ Though the specific initial conditions of equally spaced particles was 
considered it is expected to hold for more general particle distributions. For example, the case 
of Poissonianally distributed particles can be shown to give the same late-time results. From the 
physical argument given above for the deficit, it might be expected that similar behavior is also 
seen in higher dimensions. 



^even > todd) 
T^{,Zvxe.n; teven) 
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Figure 7: The late-time scaled density profile near the boundary. 



(a) Particles with v=-c 
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(b) Particles with v=+c 
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Figure 8: The late-time scaled density of (a) left-moving particles with v = — c and (b) right-moving 
particles with v = +c. 
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QB{t) 


QD(z,t) 


C(z,t) 


K(z,t) 


z < ct 


1 








1 







y/irct 




V2TTCty/l + z/ct 


z > t 


1 







i 


\J net 


\Z4TTCt 


\ZAwct 



Table 2: The late-time scaling forms for the bulk, deficit and left and right moving density profiles. 



Discussion 



In this paper we have considered the effect of an impenetrable boundary on the single-species reac- 
tion A + A — ► O. The boundary breaks the translational invariance of the system and introduces a 
spatial dependence into the reactant density profile. Two different dynamics for the reactants were 
considered: diffusive and ballistic propagation. For the case of reactants that diffuse it was found 
that an excess density is formed, whereas for ballistic annihilation a density deficit was obtained. 
It is interesting, form the view-point of surface critical phenomena to compare the influence of the 
boundary at different length and time scales on these two reaction processes. In the system with 
diffusive reactants it was seen in the RG treatment that the impenetrable/reflecting boundary 
behaves the same way regardless of scale: the propagator is not renormalized and maintains the 
form However, this is not the case for ballistic reactions. At late times the reactants tend to 
move in groups that share the same velocity. When such a group meets the reflecting boundary 
almost total internal annihilation will occur. Hence, contrary to the case of diffusive reactions, 
for ballistic annihilation the reflecting boundary becomes an effectively absorbing boundary in late 
times. 

An obvious project for further study is the effect of an impenetrable boundary on the two- 
species reaction A + B — ► O. Also, the universal properties of ballistic annihilation should be 
further explored. It is surprising to note that even in unbounded systems the problems of ballistic 
annihilation in higher dimensions, finite reaction rate and different initial velocity distributions 
have been relatively untouched. It should be possible to solve the model that was examined in 
section 2 of this paper with a finite reaction rate, using the same technique as [27\ and the same 
long-wavelength behavior is expected. More interesting would be to analyze the sensitivity of the 
evolution to different initial velocity distributions, particularly in higher dimensions. Finally, given 
the results presented in |33| it would be worth examining the effects of a rough boundary on a 
simple reaction process. 
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